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ABSTRACT 



Variability of 226 AGNs in the near-infrared J, H, and K' bands is analyzed 
and discussed. An ensemble average for measured variabilities was obtained for 
various samples of the AGNs divided by absolute _B-magnitude M B , redshift z, 
and radio strength. All the samples in the J, H, and K' bands are found to 
give significant ensemble variability, but no significant wavelength dependence 
is found. The ensemble variability in the entire sample combining the J, H, 
and K' samples is Am xs 0.22 mag, while Am 0.18 mag for the radio-quiet 
AGNs and Am 0.26 mag for radio-loud AGNs. The ensemble variability 
for the radio-quiet AGNs shows no significant M B -dependence, while showing 
positive M B -dependence for the radio-loud AGNs. In any samples the measured 
variability shows positive correlation among different passbands, with the 
correlation coefficients of tjh, thk'i and tjk> ranging from 0.6 to 0.9. For 
radio-quiet AGNs, the coefficient r HK > in a redshift range of 0.1 < z < 0.3 is 
significantly higher than r JH or r JK >. The coefficient for the radio-loud AGNs 
with 0.6 < z < 1.0 is as high as 0.95, irrespective of the passband. However, 
for the radio-quiet AGNs with z > 0.3 and radio-loud AGNs with z < 0.3, 
we cannot confirm such strong correlation among different passbands. All the 
features of near-infrared variability for the radio-quiet AGNs are consistent with 
a simple dust reverberation model of the central regions of AGNs. However, the 
features for the radio-loud AGNs are not fully explained by such a model, and a 
non-thermal variable component is suggested as a viable candidate for causing 
their large and fast variability in the near-infrared region. 



Subject headings: galaxies: active — quasars: general — galaxies: photometry 
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1. Introduction 

Active galactic nuclei (AGNs) at cosmological distance emit enormous amount of 
energy from their central region which is compact and spacially unresolved. Observations of 
AGN variability are important for the understanding of the physical mechanisms of energy 
emission by the central engine. 

Monitoring observations of many AGNs in the optical region have been made by 
various authors, and the relations among AGN optical variability, luminosity, redshift, time 
scale and so on have been derived. For example, a negative correlation between variability 
amplitude and luminosity was found, and the possibility of explaining the emission and 
variability of AGN by a sub-unit model has been discussed (e.g., Cristiani et al. 1996; Hook 
et al. 1994). The wavelength dependence of variability was discussed by various authors 
(e.g., Cristiani et al. 1997; Winkler 1997; Winkler et al. 1992). Multi-color monitoring 
observations of AGNs are useful to distinguish the dependence of valiability on wavelength 
and redshift. Optical variability was compared with ultraviolet variability (Clemente et al. 
1996), and it was concluded that the variability was larger in the shorter wavelength UV 
region than in the optical region. 

Further understanding can obviously be made by combining the near-infrared (NIR) 
data with the UV/optical data. Barvanis (1992) analyzed the optical and NIR light curves 
of Fairall 9 (Cravel, Wamsterker & Glass 1989), and explained the delay of NIR variability, 
in terms of a dust reververation model, in which thermal re-radiation comes from a hot 
dust torus illuminated by the central emission engine. Neugebauer et al. (1988) monitored 
108 PG quasars in the J, H, K, L and 10/iin bands for about 20 years. It was suggested 
that the detection rate of variability is smaller in the NIR than in the optical, if the same 
level of accuracy is required. 

We present the new data of J, H, and K' variabilities for 226 AGNs. The sample 
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was divided into groups by radio strength, absolute -B-magnitude, and redshift, and the 
ensemble variability for each group was measured. 

2. Data 

Our analysis in this paper was made by using the NIR images of AGNs obtained 
through the procedure of image reduction discussed in Paper I. Differential photometry was 
employed to measure their NIR variability. The achieved accuracy was significantly higher 
than the accuracy from an alternative method of standards-based photometry, as described 
in Paper II. Only the measurements made with more than two reference objects and having 
an accuracy higher than 0.1 mag are used in this paper. 

The AGNs in our sample were selected from various versions of the Quasars and Active 
Galactic Nuclei catalog (VV catalog, Veron-Cetty and Veron 1993, 1996, 1998). AGNs were 
selected with consideration of their use in the MAGNUM Project (Kobayashi et al. 1998a, 
1998b). The distribution of declination and right ascention for all AGNs in the sample is 
shown in Fig. 1 of Paper I. The distribution of absolute -B-magnitude and redshift is shown 
in Fig. 2 of Paper I. 

All observations were made with the 1.3m infrared telescope at the Institute of 
Space and Astronautical Science (IS AS), Japan, equipped with the NIR PICNIC camera 
(Kobayashi et al. 1994). The AGNs and reference objects were imaged in the J, H and K' 
bands with the telescope stepped in a raster pattern. Two photometric standard stars with 
different elevations were observed three times each night (for details see §2 of Paper I). 

The PICRED software, that was developed for the PICNIC camera, was used to reduce 
the images. During our observations this software was optimized for AGN and quasar 
images, to achieve fully automated reductions (for details see §3 of Paper I). 
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3. Analysis and Discussion 

3.1. ensemble variability 

Since each AGN was observed on two different nights separated by a year or more, 
the variability thus obtained does not necessarily reflect the characteristic amplitude of 
intrinsic variability. Therefore, we discuss the ensemble variability of AGNs which reflects 
the dispersion of their individual variabilities. However, the standard deviation of such data 
is not a good parameter in estimating the AGN variability, because not only the intrinsic 
AGN variability but also measurement errors broaden the distribution. 

We here introduce the ensemble variability, after excluding the contribution of 
measurement errors, as 

Am = 



N 



and its error given by 

1 



(EfAmf-Ej) (4Amfr? - 2g f ) 

M n 3 iv 2 ■ 1 ] 



2Am 

Derivation of Am and a^ m is described in Appendix [A]. In the remainder of this paper these 
quantities of Am and a^m are used to discuss the variability of the AGNs in our sample. 



3.2. The relation between the variability and AGN character 

We examine whether the variability is correlated with the parameters such as radio 
strength, rest-frame time interval of observations, absolute -B-magnitude, redshift, Seyfert 
type, and NIR colors. First, the sample was divided by each parameter into two groups, 
"a" and "6" , at a point where there appeared to be a boundary on either side of which the 
data are separated. Then, the ensemble variabilities for these two groups were compared 
with each other. 
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The AGNs in our sample are distinguished by the radio strength or the ratio of radio 6 
cm flux relative to optical V-band flux f u (6cm)/f u (V). The AGNs with f u (6cm)/f u (V) < 10 
were classified into the radio-quiet group, and others with f u (6cm)/f l/ {V) > 100 into 
the radio- loud group. The AGNs, observed twice separated by time interval At obs , were 
distinguished by the rest-frame time interval At rest = At obs /(l + z). The AGNs with 100 
days < At rcst <400 days were classified into the short-At rest group, and others with 400 
days < At rest < 800 days into the long-At rcst group. The boundary at 400 days reflects the 
period in which the AGNs were observed. Our observational runs consist of three periods 
(January 1996- April 1996, November 1996-February 1997, December 1997- April 1998). 
The AGNs observed in the first and third periods are mainly of long interval, and those 
in the second and third periods are of short interval. The AGNs with M B < —23.5 were 
classified into the bright group, and those with M B > —23.5 into the faint group. The 
AGNs with z < 0.3 are classified into the low- z group, and those with z > 0.3 into the 
high-z group. The AGNs with Seyfert 1, 1.2, and 1.5 were classified into the early- type 
Seyfert group, and others with 1.8, 1.9, and 2 into the late-type Seyfert group. The AGNs 
with J — H < 0.8 and those with H — K' < 0.8 were classified into the blue group, and 
others with J — H > 0.8 and H — K' > 0.8 into the red group. 

The ensemble variabilities for the "a" and "6" groups are derived for each of radio 
strength, rest-frame time interval of observations, absolute -B-magnitude, redshift, Seyfert 
type, and NIR colors. Table 1 shows the ratio Am("o")/Am(T) for each of the above 
parameters. The last column of this table represents the average ratio taken over the J, H, 
and K' bands. The average ratio between the radio-loud and radio-quiet groups is 1.46, 
which is the largest. The average ratio for the long- and short-At rcst groups is 1.24, while 
it is 1.20 for the bright- and faint-M# groups. The average ratios for other quantities are 
much closer to unity. 
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Top panel of Fig. 1 shows that the ensemble variability in the J, H, or K' band is 
Am ~ 0.2, estimated for the entire sample, while Am 0.25 — 0.3 and 0.18 — 0.2 for the 
radio-quiet and radio-loud groups, respectively. We see no significant wavelength-dependence 
of Am. 

We divide the radio-quiet or radio-loud group furthermore into two subgroups of short 
and long At rest , and estimate the ensemble variability for respective subgroups. Bottom 
panels of Fig. 1 show that the ensemble variability in the J, H, or K' band is Am 0.2 
(long At rest ) and 0.15 (short At rcst ) for the radio-quiet sample, while Am ~ 0.28 — 32 (long 
At r est) and 0.22 — 0.28 (short At rcst ) for the radio-loud sample. 

We furthermore divide the short- or long-At rest subgroup by Mb or z. In this way, the 
ensemble variabilities of the radio-quiet AGNs in respective subgroups are shown in Fig. 2. 
The similar results for the radio-loud AGNs are shown in Fig. 3. 

3.2.1. statistical test and estimation of the A dependence of the ensemble variability 

A statistical test on the wavelength dependence of Am^ was done by applying the 
X 2 method to the result. Since such dependence was not found to be significant, it is 
reasonable to adopt a two-parameter function of Am\(ai,a2) = a\e a2X and search for the 
solution near a 2 ~ in minimizing x 2 — J2\=j,H,K>(Am\ ~ Am^ai, a 2 )) 2 / cr 2 Xmx . Table 2 
shows the optimized values of a 1 = 0.205 and a 2 = 0.029 for the entire sample, leading to 
Am K ,/Amj = 1.03, otherwise Am K , / Amj = 0.83 - 1.26 (95% C.L.) and 0.79 - 1.32 (99% 
C.L.). For the radio-quiet sample, the optimized values of a-y = 0.205 and a 2 = —0.065 give 
Am k' I Am j = 0.94. A similar result holds if the radio-quiet sample is further divided by 
At rcst . For other groups, the C.L. range becomes wider, so that the ability of rejecting the 
hypothesis of no wavelength-dependence by the test remarkably decreases. 
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3.2.2. statistical dependence test of the ensemble variability on characteristic parameters 

In this section, using the ensemble variability Am averaged over the J, H, and 
K' bands, we examin whether Am depends on radio strength, At rest , M B , and z. 
In this case, with a simple one-parameter function of Am(ai) = a±, it is reduced to 
X 2 = J2j(Amj — Am(ai)) 2 / a\ mj . Table 3 shows the optimized value of a\ together with P 
which represents the reliability of rejecting the hypothesis that Am does not depend on the 
parameter in question. 

It is understood from this table that the statistical equivalence of Am between the 
radio-quiet and radio-loud samples is rejected by a level of P > 99.9%. For the radio-quiet 
sample further divided by At rcst , the statistical equivalence of Am between the short and 
long At rest is also rejected by a level of P > 99.9%. For the radio-loud sample, however, 
such statistical equivalence is rejected only by P = 71.1%. 

For the radio-quiet sample with short At rest , the test for the M B -dependence gives 
P = 75.0% by which it is difficult to conclude with certainty that Am depends on M B . The 
same test for z-dependence gives P = 95.5%, indicating a rather strong z-dependence of 
Am. For the radio- loud sample with short At rcst , both M^-dependence and ^-dependence 
are highly significant. We note that such strong M^-dependence is in clear contrast with 
the result for the radio-quiet sample, and such strong z-dependence is similar to the result 
for the radio-quiet sample. 

3.2.3. discussion of the parameter dependence of the ensemble variability 

Figure 4 shows the relations among various parameters for the radio-quiet sample. In 
the left column, from top to bottom, are shown the values of Am, At rcst , and Mb, estimated 
in four z-bins between z = and 1. Similarly, in the right column, from top to bottom, 
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are shown the values of Am, At rest , and z, estimated in four M B -bins between M B = —32 
and —18. In this figure, open and filled squares represent the samples with short and long 
At rest , respectively. 

We see from the left column that Am, At rest and M B decrease with increasing z, except 
for the case of Am in the lowest-z bin. These trends with z are equivalently converted to 
the trends with Mb, as shown in the right column, by using the monotonical z versus Mb 
relation. We caution that the z-dependence of At rest may only be apparent, arising from 
the cosmo logical time delay At rcst = At Q b s / (1 + z) applied to our sample which has a rather 
limited range of At b s - Consequently, any trends with At res t may also be apparent. 

The ensemble NIR variability of radio-quiet AGNs is as small as Am < 0.25, showing 
little Ms-dependence in this work. We note that faint AGNs, mostly at low z, are 
contaminated by a host galaxy component, which is indicated by our mult i- aperture color 
analysis (Paper I; see also Kotilainen & Ward 1994). Since the host galaxy component is 
stellar and does not vary on a time scale of years, such contamination has the systematic 
effect of weakening the AGN variability. Therefore, Am, after correction for this effect, 
would still have little or negative correlation with absolute _B-luminosity, depending on the 
degree of contamination within chosen aperture size. 

The above M B -dependence of Am is expected from the model of dust reverberation in 
which brighter AGNs have a larger dust torus. That is, a variation in the UV/optical light 
emitted from the central engine is absorbed in more extended region of dust from which 
the NIR radiation is emitted and the spread in arrival times of the NIR variation from 
this extended region produces a variation with smaller amplitude. Therefore, regardless of 
the real emission mechanism of central source, the NIR variability would show only small 
correlation with M B , as observed. 

Next, we consider the radio-loud AGNs. Figure 5 shows their relations among Am, 
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At rest , M B , and z, in a similar way as in Fig. 4. We divide the radio-loud sample into the 
short and long At rest . However, because of the lack of enough data, the values of Am and 
At res t are not estimated as a function of z or Mb, for the case of long At rest . 

We notice that Am for the radio-loud AGNs, in the case of short At rest , strongly 
increases with increasing z or with increasing absolute ^-luminosity, opposite to that for 
the radio-quiet AGNs. If we assume that such positive correlation is only apparent, wishing 
to explain it in terms of the effect of contamination of host galaxy component, we have to 
invoke an extremely different contribution from the host galaxies between radio-loud and 
radio-quiet AGNs, which is difficult to justify. Therefore, it is more reasonable to conclude 
that the different Ms-dependence of Am reflects the different emission and variability 
mechanisms between radio-loud and radio-quiet AGNs. 

3.3. The correlation of variability in different passbands 

In this section the correlation among the variabilities in the J, H, and K' bands is 
discussed using only the data of estimated variabilities with more than two reference stars 
and an accuracy better than 0.1 mag. The results for the radio-quiet and radio-loud AGNs 
are shown in Figs. 6 and 7, respectively, where Am H is plotted against Amj, Am' K against 
Amn, and Am' K against Amj. Open and filled symbols correspond to long and short 
At rest , respectively. 

The data distribute, more or less, along the diagonal running from lower left to upper 
right through the origin, which indicates that AGNs, becoming brighter in one band, 
become brighter in the other band, and vice versa. It is seen from each panel that more 
data are plotted in the lower left region than in the upper right region. However, this is not 
real because AGNs, becoming fainter, are likely to be either undetected or rejected by our 
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accuracy requirement. 

The correlation coefficient of variabilities in two different bands, say J and H, is 
defined as 

E(Amj,i - Am,j)(Am Hti - Am H ) , . 

fjH — —, = , 1<jJ 

VElAmj,, - A mj ) 2 E(A m//i! - Am H ) 2 

where Am x (A = J or H) is the unweighted average of variability data in each band. Table 
4 tabulates the values of r JH = +0.74, r H K' = +0.81, and r JK / = +0.71 for the entire 
sample, for which the 68.3% confident interval is about 0.1 or less. The coefficient is higher 
than +0.59 for the radio-quiet sample, and even higher than +0.8 for the radio-loud sample. 

Figure 8 shows the z-dependence of rjn, rHK', and tjk 1 - For the radio-quiet sample, 
all these coefficients at z < 0.1 are equally high. While r JH keeps a high value irrespective 
of z, the coefficients r JK i and r HK i are getting smaller for higher z. This trend may have 
occurred from an underestimation of the correlation, because of lower statistical accuracy in 
the high-2 sample of smaller size, and because of larger errors in Am for the high- z sample 
consisting mainly of faint AGNs. For the radio-loud sample, all the coefficients keep a high 
value irrespective of z and their accuracy becomes higher for higher z, in sharp contrast 
with the radio-quiet sample. 

The equivalence of r JH , r HK /, and r JK i was tested against the estimated true value of 
coefficient r trnc . Table 5 tabulates the estimation of r truc together with P which represents 
the reliability of rejecting the hypothesis that rjn, thk', and tjk> are equivalent to each 
other. For the radio-quiet sample, their equivalence at z — 0.1 — 0.3 is significantly rejected, 
while not definitely so at higher z. For the radio-loud sample, the equivalence at z > 0.3 is 
rejected with negligible reliability, in other words, it is very probable that r JH , r HK >, and 
tjk> are the same at z > 0.3. 

It is important to note that rjn, thk 1 , and tjk> for the radio-quiet sample are equally 
high at z < 0.1, and thk' is higher than tjh and tjk> at z — 0.1 —0.3. In dust reverberation 
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model, the NIR flux of AGNs is mainly emitted from hot dust which is heated up to the 
evaporation temperature T evap xs 1500K (e.g., Kobayashi et al. 1993). The equivalence 
of tjh, Thk'i and tjk> at z < 0.1 can be explained by the black body radiation with 
constant T evap ~ 1500K dominating on the longer-wavelength region of the 1/xm minimum 
seen in the rest-frame SED of AGNs and quasars. On the other hand, the main variable 
component dominating on the shorter- wavelength region of the 1/zm minimum is considered 
to be a power-law component, and its variability is not necessarily synchronized with the 
NIR variability. Figure 9 shows this non-synchronization between UV/optical and NIR 
variabilities on either sides of 1/im in rest frame. Thereby, at z — 0.1 — 0.3, the correlation 
of variabilities between the J band and longer wavelengths becomes weak, because the 1/xm 
minimum in the rest frame moves to the J band and the power-law variable component 
affects the flux there. At z > 0.3, the 1/zm minimum moves in between the H and K' 
band. In such case, tjh is kept high, determined mainly by the the power-law variable 
component, while thk 1 and tjk> become low, determined by both the power-law and 
black-body components. 

For the radio-loud sample, however, the obseved trends of r JH , r HK >, and r JK i are very 
different from those for the radio-quiet sample and may be explained by a mechanism of 
variability other than dust reverberation. 



3.4. The time scale of AGN variability 

In this section we consider the relation between the time scale and other parameters 
that characterize the variability of AGNs. For pedagogical purpose, we introduce two 
representative functions to be fitted to our NIR data: 

Am(A,p) = A(At rcst f , (4) 
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and 

Am(B,r) = - exp(-At rcst /r)) . (5) 

These functions correspond to a divergent increase of variability and an asymptotic increase 
with increasing At rest , respectively 

Tables 6a and 6b tabulate the fitted values of the parameters and their 95% and 99% 
confidence intervals, applied to the AGNs in the entire sample as well as the radio-quiet 
and radio-loud samples. The curves with such fitted values are shown in Fig. 10. We note 
that p and r are searched for in their postive values to avoid division by at At rcst = 
in equations 4 and 5. If their fitted values or the lower limits of their confidence intervals 
became negative, they were given the value of zero. 

All the fitted curves are accepted independent of how the data are divided into 
subgroups. Thus, it is difficult to prefer the function in equation 4 to the other in equation 
5. However, it is generally suggested from Tables 6a and 6b that p and r for the radio-loud 
sample are smaller than those for the radio-quiet sample, which indicates that the time 
scale of NIR variability of radio-loud AGNs is shorter than that for the radio-quiet AGNs. 

Figure 1 1 shows the fitted values of p and r in comparison with those obtained by 
Cristiani et al. (1996) from the optical variability data in the literature. The vertical 
errorbar is for the 68.3% confidence interval, and the horizontal errorbar for the covered 
range of wavelengths. It is interesting to see that the time scale of optical variability by 
Cristiani et al. (1996) agrees, within the range of uncertainties, with the time scale of 
NIR variability for the radio-quiet sample. On the contrary, this is not the case for the 
radio-loud sample, where the time scale of optical variability is significantly longer than the 
time scale of NIR variability. 

From the viewpoint of dust reverberation, it is natural that the NIR variability occurs 
on a rather elongated time scale, because the UV/optical variability is re-emitted in the 
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NIR from a more extended region of dust producing a spread of arrival times. Consequently, 
whether the NIR variability occurs as fast as or faster than the UV/optical variability gives 
a clue which limits or rejects a mechanism of variability based upon the dust reverberation 
model. In this regard, it is difficult to explain by the dust reverberation model how the NIR 
variability of the radio-loud AGNs occurs faster than the optical variability, as seen from 
the comparison in Fig. 11. 



3.5. The relation between near infrared variarility and radio activity 

We have highlighted the difference of NIR variability between the radio-quiet and 
radiou-loud AGNs. Obtained features for these respective samples are summarized as 
follows: 

(1) The amplitude of ensemble NIR variability Am in the radio-quiet sample is smaller 
than in the radio-loud sample. 

(2) The wavelength-dependence of Am in the NIR region is not found in both the radio-quiet 
and radio-loud samples. 

(3) No correlation is found between Am and Mb in the radio-quiet sample, but a negative 
correlation is suggested if corrected for the effect of possible contamination of a host galaxy 
component. On the other hand, a positive correlation is found between Am and Mb in the 
radio- loud sample. 

(4) The correlation coeficient thk between the H and K variabilities at z — 0.1 — 0.3 is 
significantly higher than tjh and tjk 1 in the radio-quiet sample. On the other hand, the 
coeficients rjn, tjk> and tuk 1 at z > 0.3 have a high value of 0.9 — 0.95 in the radio-loud 
sample. 
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(5) The time scale of NIR variability in the radio-quiet sample is longer than in the 
radio-loud sample. 

The features of 2, 3, and 4 for the radio-quiet AGNs are consistent with those expected 
from a mechanism of variability by dust reverberation. Furthermore, the feature of 5 for 
the radio-quiet AGNs, if combined with the time scale of optical variability by Cristiani et 
al. (1996), is also consistent with dust reverberation. 

However, the features of 2 to 5 for the radio-loud AGNs are not explained consistently 
by dust reverberation. A non-thermal variable component, as a substitute of thermal 
radiation from hot dust, is worth considering. For example, Sanders et al. (1989) proposed 
such a non-thermal component which would vary more strongly and faster in the compact 
region than the thermal radiation from an extended hot dust torus. 

We proceed with the working hypothesis such that dust reverberation is responsible for 
the emission and variability of radio-quiet AGNs, while a non-thermal variable component 
is responsible for the emission and variability of radio-loud AGNs. In the next subsections, 
we check the plausibility of this hypothesis using additional data from the literature. 



3.5.1. The radio strength of the dust reverberation sample 

It is known that the time delay of NIR variability lagged behind the the UV/optical 
variability is a key prediction of the dust reverberation model. Therefore, measurements 
of UV/optical and NIR light curves, based on multi- wavelength monitoring observations 
of AGNs, immediately reject or justify the application of the dust reverberation model to 
individual AGNs. 

The AGNs with measured time delay, to which the dust reverberation model 
is successfully applied, are taken from the literature, and their data of Mb, and 
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f„(6cm)/ fviV) are summarized in Table 7. It is evident that these AGNs are classified as 
radio-quiet with f v {<6cm)/ f v (V) < 10. The sample of these AGNs is heterogeneous and 
is far from complete. It is possible that the sample is biased in favor of nearby and/or 
radio-quiet AGNs. 

The first systematic monitoring observations of many Seyfert 1 AGNs and quasars in 
the V and K bands were performed by Nelson (1996a). Table 8 shows the statistics based 
on Nelson's sample consisting of 51 program AGNs. The first row represents the numbers 
of radio-quiet and radio-loud AGNs in his entire sample of 51 AGNs, the second row for 
the subsample of 33 AGNs which were found to vary in both the V and K bands, and the 
third row for the subsample of 6 AGNs for which the time delay of the -fT-band variability 
relative to the y-band variability were measured. 

We see that the AGNs with measured time delay, which are best explained by the 
simple dust reverberation model, are always classified as radio-quiet. In other words, 
radio-quiet AGNs are potential targets for mult i- wavelength monitoring from which the 
time delay between the NIR and UV/optical variabilities can certainly be measured. 

3.5.2. The relation between variability and flatness of radio SED 

The large and fast NIR variability of radio-loud AGNs found in this paper is not 
explained by dust reverberation model. In order to understand what causes such variability, 
we further classify the radio-loud AGNs with respect to their spectral feature in the radio 
region, being either flat (a > —0.5) or steep (a < —0.5), if fitted to a power-law form of 
f v oc v a (Sanders et al. 1989). 

The values of power index a in our sample were determined using the data of radio 
fluxes at 6 cm and 11 cm taken from the VV catalog. Figure 12 shows that a is not 
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correlated with M B or z. Figure 13 shows that the ensemble variability in the J, H, 
and K' bands is Am « 0.35 — 0.45 for the radio-loud and flat-spectrum AGNs, which is 
systematically higher than Am R3 0.2 for the radio-loud and steep-spectrum AGNs. It is 
therefore the flat component that actually brings about the large and fast NIR variability. 

The radio spectral index, a, is known to be well correlated with the radio-loud AGN 
morphology; while the steep spectrum is associated with extended lobe-dominant sources, 
the flat spectrum is associated with core-dominant variable sources such as OVVs, highly 
polarized QSOs, and BL lacs. Thereby, it is reasonable to conclude that a non-thermal 
variable component, as exemplified by non-thermal emissions from such objects (e.g., 
Robson et al. 1993; Bloom et al. 1994), is responsible for occurrence of features of the NIR 
variability found for the radio-loud AGNs in this paper. 



4. Summary 

We presented comprehensive study of NIR variability of 226 AGNs based on multiple 
observations in the J, H, and K' bands. Our sample consists mainly of Seyfert 1 AGNs and 
QSOs. About a quarter of objects in each category are radio loud. The AGNs in the entire 
sample have redshifts spanning a range from z — to 1, and the absolute -B-magnitudes 
from M B = -29 to -18. 

Based on the method of differential photometry, we find that the ensemble NIR 
variability for the entire sample of AGNs is typically Am 0.2 mag. When the sample 
is divided by radio strength, the variability for the radio-quiet sample is systematically 
smaller than that for the radio-loud sample. No clear wavelength-dependence of Am in the 
NIR region is found for either the radio-quiet or radio-loud sample, in sharp contrast with 
the UV/optical result in the literature. 
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We examined the dependence of Am on various quantities such as radio strength, M B , 
and z, with special attention as to whether their ensemble variability would support or 
reject the simple-minded dust reverberation model for AGNs. 

The radio-quiet AGNs show no significant correlation between Am and M B) although 
negative correlation is suggested if corrected for the effect of possible contamination by 
a host galaxy component. On the other hand, the radio-loud AGNs show a positive 
correlation between Am and M B . 

The radio-quiet AGNs give a significantly higher correlation coefficient r HK between 
the H and K variabilities at z — 0.1 — 0.3, when comapred to tjh and tjk>- On the other 
hand, the radio-loud AGNs give a high value of 0.9 — 0.95 to all the coefficients r JH , r JK i 
and thk 1 at z > 0.3. 

Time development of ensemble variability is examined using heuristic functions. The 
time scale of NIR variability for the radio-quiet AGNs in this paper is not shorter than 
the time scale of UV/optical variability given in the literature. However, the time scale of 
NIR variability of the radio-loud AGNs is significantly shorter than the time scale of their 
UV/optical variability. 

All the features of NIR variability for the radio-quiet AGNs are consistent with those 
expected from the dust reverberation model. However, it is difficult for this simple-minded 
model to explain the features of the radio-loud AGNs, and a non-thermal variable 
component is suggested as a viable candidate for causing the large and fast NIR variability 
of the radio-loud AGNs. 
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A. estimation of the ensemble variability 

This section describes how we estimate variability and its error in this work using the 
data for which the sample number and accuracy are both limited. The data obtained in 
this work are given in the format of (Ami, ax), (Am2,f2) ••• {Am^ -,<Jn) for the sample of 
AGNi, AGN 2 ... AGNjv- Here, Awij and <7j represent the variability and its error of the z-th 
object. Hereafter, the real variability of the i-th object is described as Vi, and the difference 
between t>i and Ami is given by e*. Then, can be regarded as a random variable described 
by the normal distribution having zero average and the standard deviation <7j, that is, 

Ami = Vi + e u < e { >= 0, < e\ >= cr h (Al) 

where <> represents the expectation. The real variability of the i-th object V{ is assumed 
to be a random variable that follows the normal distribution with zero average and the 
standard deviation v independent of i. 

In the ideal case where the data have zero error and the sample size is infinite, the 
dispersion of variability is given by 

N „,2 



lim ^-p- (A2) 

Hereafter, the aim is to estimate the value of equation A2 and the error from the existing 
data. Consider 

v] = Am\ - 2vi€i - el (A3) 
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Then, the expectation of the second term in the right-hand side vanishes because V{ and 
are independent of each other. Therefore, the expectation of v 2 becomes 

<v 2 >=Am 2 -a 2 . (A4) 

Substitution of the observed data into equation [A^ gives 



lim 

N^oo 



N 



T N v 2 V N Am 2 -V N a 2 — 

UJ\ „ u ^ u °> = vl (A5) 

where the first ~ stems from the finite sample number and the second ~ stems from the 
error of the data. The overline of v 2 indicates that the value is the expectation based on 
the exisiting data. 



The error of equation |A5| is estimated as follows: First, we focus on the error caused 
by the deviation of the ideal expectation from that of the actual data with finite sample 
number. It should be noticed that Y^f {vi/vo) 2 is the random variable described by the x 2 
distribution with N degrees of freedom. Since Vi/vq (i = 1, 2...N) is a random variable that 
follows the normal distribution with zero average and the standard deviation of unity, the 
expectation X)f (ui/i>o) 2 is equal to iV with the dispersion of 2N. Therefore, the error of 



equation arising from the finite sample number is estimated as ±^2vq/N. 

We next focus on the error caused by the deviation of the ideal expectation from that 
of the actual data with non-zero error. Assuming that V{ and e» are independent of each 
other, the difference between v 2 and < v 2 > is calculated as 

<{v 2 -<v 2 >) 2 > = <(Am 2 -2e l v i -e 2 -Am 2 + a 2 ) 2 > 

= < Ae 2 v 2 + e\ + af + 4efu 4 - 4e^cr 2 - 2e 2 a 2 > 

= 4a 2 <v 2 > + <et>-af, (A6) 

where equation pO is used. Using equation lAl and < ef >= 3af, equation RE furthermore 



reduces to 



< (vf- < vf >) 2 > = Aaf < vf > +2af 
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= 4a 2 Am^-2af. (A7) 
Therefore, the error of equation |A5| arising from the non-zero error of the variability data is 



estimated as ±vEi (4Am?of — 2crf)/N 2 . Using equation |X5| and the above error estimates 



the expectation value of Vq and its error are estimated as 



- Ef Am?-£jV , 4 M Ef (4Am? ( x| - 2*?) 
w o = r T ± \ TT + 



V N N 2 



iv ± V + • (A ^ 

Finally, using the usual error propagation, the expectation value of vq and its error are 
estimated as 



N 
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Fig. 1. — The ensemble variability of AGNs in the J, H and X' bands. Shown are the results 
for all AGNs, and radio-quiet and loud AGNs (upper panel), for radio-quiet AGNs divided 
by At rcst (lower left panel), and for radio- loud AGNs divided by At rcst (lower right panel). 
Only the data with more than two reference objects and with accuracy higher than O.lmag 
are plotted. 

Fig. 2. — The ensemble variability of radio-quiet AGNs in the J, H and K' bands. Shown 
are the results for short- At rcst AGNs divided by M B (upper left panel) or by z (lower left 
panel), and those for long-At rcst AGNs divided by Mb (upper right panel) or by z (lower 
right panel). 

Fig. 3. — The ensemble variability of radio-loud AGNs in the J, H and K' bands. Shown 
are the results for short- At rcst AGNs divided by Mb (upper left panel) or by z (lower left 
panel), and those for long-At rcst AGNs divided by Mb (upper right panel) or by z (lower 
right panel). 

Fig. 4. — The relation among characteristic parameters for radio-quiet AGNs. Shown are 
Am, At rest , and Mb plotted against z (left column), and Am, At rcst , and z plotted against 
Mb (right column). 

Fig. 5. — The relation among characteristic parameters for radio-loud AGNs. Shown are 
Am, At rest , and M B plotted against z (left column), the Am, At TCSt , and z plotted against 
Mb (right column). 

Fig. 6. — The relation between variabilities in two NIR bands for radio-quiet AGNs divided 
by At rest . 

Fig. 7. — The relation between variabilities in two NIR bands for radio-loud AGNs divided 
by At rest . 
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Fig. 8. — The coefficient of correlation between AGN variabilities in the \ and Xj bands 
as a function of redshift z. Shown are the results for radio-quiet AGNs (upper panel), and 
for radio-loud AGNs (lower panel). 

Fig. 9. — The ensemble variability of AGNs as a function of rest-frame frequency. Results 
in the NIR are from this paper, and those in the UV/optical are from the literature. Filled 
small symbols are for 3 month-data (triangles, Cimatti et al. 1993; filled circles, Trevese et 
al. 1994; squares, De Clemente et al. 1996; pentagons, Hook et al. 1994; hexagons, Cristiani 
et al. 1990). The open small symbols are the same as filled ones, but for 2yr-data. 

Fig. 10. — The ensemble variability of AGNs as a function of rest-frame time interval between 
two observations. Shown are the results for all AGNs (upper panel) and for radio-quiet 
and loud AGNs (lower panel). The solid and dotted lines represent the best-fit curves of 
Am(A,p) = A(At TCSt ) p and Am(B,r) = B(l — exp(— At rest /r)), respectively. 

Fig. 11. — The wavelength-dependence of fitted parameters to the ensemble variability of 
AGNs, such as the power index p in Am(A,p) = A(At rcst ) p ( Upper panel), and the timescale 
r in Am(B,r) = B(l — exp(— At rest /r)) (Lower panel). The vertical errorbars represent 
the 68.3% confidence interval, and horizontal errorbars represent the range of wavelength 
coverage in the analysis. Results in the NIR are from this paper, and those in the optical 
are from Cristiani et al. (1996). 

Fig. 12. — The distribution of radio-loud AGNs in the a versus Mb diagram (upper panel) 
and in the a versus z diagram (lower panel). The boundary between the steep and flat 
spectra is placed at a = —0.5. Open and filled circles correspond to short- and long-At rest 
samples, respectively. 

Fig. 13. — The ensemble variability of radio-loud AGNs in the J, H and K' bands. Shown 
are the results for those with flat and steep radio spectra. 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 7 
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Table 1. Ratio of ensemble variabilities in two groups of AGNs 



Am("a")/Am(T) 



Tested property 


J 


H 


K' 


Average 


Radio strength 


1.37 


1.37 


1.65 


1.46 




1.22 


1.37 


1.14 


1.24 


M B 


1.29 


1.18 


1.14 


1.20 


Redshift 


1.05 


0.95 


1.16 


1.05 


Seyfert type 


1.06 


1.04 


0.92 


1.01 


J-H 


1.02 


1.00 


1.01 


1.01 


H-K' 


1.25 


1.07 


0.99 


1.10 



Note. — The entire sample is divided by each parameter 
into the "a" and "6" groups, such as (/„(6cm)//„(V) > 100, 
M6cm)/.U(V) < 10), (400 < At rcst < 800, 100 < At rost < 
400), (M B < -23.5, M B > -23.5), (0 < z < 0.3, z > 0.3), 
(Seyfert 1-1.5, Seyfert 1.8-2), (J-H < 0.8, J-H > 0.8), 
and (H-K' < 0.8, H - K' > 0.8). 
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Table 2. Test for wavelength dependence of the ensemble variability 



Confidence interval of a-i 
Sample a x a 2 95% 99% 



All AGNs 0.205 0.029 -0.192 - 0.242 -0.248 - 0.295 

Radio-quiet AGNs 0.205 -0.065 -0.333 - 0.193 -0.403 - 0.258 

Short Ai rost 0.139 -0.044 -0.575 - 0.452 -0.741 - 0.593 

Long Ai rost 0.242 -0.079 -0.404 - 0.234 -0.494 - 0.316 

Radio-loud AGNs 0.213 0.136 -0.342 - 0.582 -0.470 - 0.693 

Short Ai rcst 0.190 0.169 -0.398 - 0.703 -0.556 - 0.841 

Long Ai rost 0.251 0.111 -0.952 - 0.992 -1.313 - 1.228 



Note. - ■ The test was done by adopting a two-parameter function of Am\ = 
a\ exp(a2A). 
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Table 3. Test for radio strength, At rest , M B , and z dependence of the ensemble variability 



Sample Tested property a\ <r ai \ 2 P{%) 



All AGNs Radio strength 0.196 0.007 15.767 >99.9 

Radio-quiet AGNs Ai rcst 0.167 0.008 27.797 >99.9 

Short At rcst M B 0.126 0.011 2.757 75.0 

z 0.115 0.010 6.179 95.5 

Long Ai rost M B 0.212 0.012 1.167 44.0 

z 0.196 0.011 14.692 >99.9 

Radio-loud AGNs At rcst 0.262 0.019 1.125 71.1 

Short At rcst M B 0.131 0.016 35.898 >99.9 

z 0.117 0.013 39.503 >99.9 



Note. — P represents the reliability of rejecting the hypothesis that AmA = a\ 
does not depend on the tested property. 
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Table 4. Correlation coefficient of variabilities in the A, and Xj bands 



Sample 


Xi, Xj 


n 




Confidence interval (68.3%) 


All AGNs 


J, H 


127 


0.74 


0.70 - 0.78 




H, K' 


111 


0.81 


0.77 - 0.84 




J, K' 


118 


0.71 


0.66 - 0.76 


Radio-quiet AGNs 


J, H 


91 


0.65 


0.59 - 0.71 




H, K' 


84 


0.72 


0.66 - 0.77 




J, K' 


85 


0.65 


0.59 - 0.71 


Radio-loud AGNs 


J, H 


31 


0.88 


0.83 - 0.91 




H, K' 


24 


0.91 


0.87 - 0.94 




J, K' 


30 


0.86 


0.80 - 0.90 



Note. - 



n represents the number of AGNs in a sample. 
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Table 5. Equivalence test for the correlation coefficients tjh, thk 1 , and tjk> 



Sample 


rjH 




rjK> 


''"true 


P(%) 


Radio-quiet AGNs 












z - 0.0-0.1 


0.69 (36) 


0.77 (38) 


0.77 (38) 


0.75 


31.5 


z = 0.1 -0.3 


0.63 (29) 


0.89 (23) 


0.57 (22) 


0.72 1 


96.8 


z = 0.3 -0.6 


0.65 (19) 


-0.04 (17) 


0.50 (19) 


0.42 


58.1 


z = 0.6-1.0 


0.52 (7) 


-0.51 (6) 


-0.11 (6) 


0.03 


68.3 


Radio-loud AGNs 












z = 0.1 -0.3 


0.74 (6) 


0.76 (5) 


0.04 (5) 


0.60 


46.0 


z = 0.3 -0.6 


0.70 (10) 


0.67 (8) 


0.72 (11) 


0.70 


1.17 


z- 0.6 -1.0 


0.94 (12) 


0.97 (9) 


0.95 (12) 


0.95 


12.9 



Note. The number in the parensatheses is the number of AGNs in a 

sample. r trU c is the estimation of true correlation coefficient. P is the reliability 
of rejecting the hypothesis that rjn, Thk', and tjk 1 are equivalent to each 
other. 

1: The value of r truo = 0.72 is not adequate since the equivalence is rejected. 
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Table 6a. Fitted values of the parameters in a function of Am = A(At TCSt ) p 









Confidence interval of p 


Sample 


A 


P 


95% 


99% 


All AGNs 


0.044 


0.26 


- 0.53 


- 0.60 


Radio-quiet AGNs 


0.0036 


0.63 


0.30 - 1.00 


0.22 - 1.09 


Radio-loud AGNs 


0.077 


0.21 


- 0.79 


- 0.91 



Note. — The cases with no lower limit to the confidence interval are those 
for which the lower limit becomes negative. 
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Table 6b. Fitted values of the parameters in a function of Am = B(l — exp(— At rcst /r)) 









Confidence interval of r 


Sample 


B 


r 


95% 


99% 


All AGNs 


0.24 


180 


50 - 370 


- 440 


Radio-quiet AGNs 


0.30 


510 


200 - 


170 - 


Radio-loud AGNs 


0.29 


140 


- 820 


- 2040 



Note. - The cases with no lower limit to the confidence interval are 
those for which the lower limit becomes negative. The cases with no upper 
limit are those for which linear relation of Am oc Ai rost is accepted from 
the x 2 test. 
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Table 7. AGNs in the literature which support a dust reverberation model. 



Object Ref z M B f v (6cm)/f v (V) 



GQ COM 1 0.165 -24.4 0.0 

Fairall 9 2 0.046 -23.0 0.0 

NGC 3783 3 0.009 -19.7 0.7 

MARK 744 4 0.010 -19.3 0.0 

NGC 4151 5 0.003 -18.7 2.3 

NGC 1566 6 0.004 -18.0 3.7 



Note. — The data of z, Mb, /;,(6cm), and f v (V) are taken 
from the VV catalog. All AGNs listed here are regarded as 
radio quiet satisfying our criterion of /„ (6cm) < 10. 
References: (1) Sitko, M. L. et al. 1993, (2) Clavel, J., 
Wamstecker, W. and Glass, I. S. 1989, (3) Glass, I. S. 1992, (4) 
Nelson, B. O. 1996b, (5) Oknyanskii, V. L. 1993, (6) Baribaud, 
T. et al. 1992. 
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Table 8. Statistics of AGNs based on Nelson's V and K band monitoring observations. 



Sample 


Number of AGNs in ; 


i sample 


Radio qiet 


Radio loud 


Ambiguous 


All AGNs 


43 


5 


3 


AGNs with detected variability 1 


29 


3 


1 


AGNs with time delay obtained 2 


6 









Note. — The result of variability and time delay is taken from Nelson (1996a). 
1: AGNs for which the variability is detected in cither V or K band or in both. 
2: AGNs for which the variability is detected in both V and K bands and the time 
delay of light variation in the K band relative to the V band is measured. 



